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APPLICATION OF DIGITAL TEMUIN DATA TO QUANTIFY AND 
REDUCE THE TOPOGRAPHIC EFFECT ON LANDSAT DATA 

I. Intfodttction 

The topographic effect » defined aa the vafiatioii in ndianccs from incl i ned ntifnees compered 
to radiance from a hoiuonta! nirface as a function of the orientation of the nufaccs relative to the 
light source and sensor position (Holben and Justice i 9h0). On Landsat unagm of rugged terrain, 
this effect is muiifested b>' the visual impression of relief (figure 1 1. Holhen and Justice ( I *179) 
measured the topographic effect on remotely sensed data and shotted the effect to be most ex* 
treme at low solar elevations and greatest for slopes in the principal plane of the sun. They also 
shosved by a Landsat simubtion study that the t«jpographic effect can produce a considerable 
vahatioii in radbnees assoebted with a given cover type and may bad to poor cover-classifkation 
results. Sadosrski and Malib ( 1977) demonstrated that reflectances vary as a function of slope 
and aspect and that such terrain variations complicate the tasL of disniminating woodland cate- 
gories with remotely sensed data. Several other studies have deni^.istrated the need to consider 
lonographic variations when undertaking Landsat cover ebssification for areas of mountainous 
terrain (Hoffe*' and Staff 1 9‘’5. Anuta 1 97b. Milbr t7 a/. 1 978. Strahbr er a/. 1978. Hoffer t-r a/. 

I 97Q. Strahbr ef al 1 97q. Williams and Miller 1 979). 

The above studies incorporated elevation data with Landsat data to improve cover classifi- 
cation accuracies. Removal or reduction of the topographic effect before classification will reduc'e 
'he variation associated with the radiance for a given cover type and thereby increase the likeli- 
hood of ebss separability. Sadowski and Maliia (1977) showed that changes in reflectance resulting 
from different topographic location resulted in spectral overlap of forested and nonforested sites. 
Williams ei al. ( 1979) demonstrated that the utility of Landsat data for mapping gypsy moth 
defoliation of forest canopies in the mountains of central Pennsylvania was severely restricted by 
topographic variations. This restriction was most pronounced in separating levels of defoliation on 
slopes of different orientations. 
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The ol^thre of this study is to quantify and reduce she topogiapluc effect found in L a ndsat 
data for an area adjacent to that examined hy Williams et W. ( I479> in central ^ennsyhrania. The 
next section of the paper describes the field area selected for this study and the methods and data 
base used in the analysis. The third section discimes four techniques for reducint the topopaphic 
effect. The fourth section examines the strength of topofraphic effect in the data and the effective* 
ness of the four techniques for reducinf the effect. The fifth section summarizes the results of the 
analysis and discusses the implications of the results. 

Oescnptioi of study area, data base, and methods 

The study area is located in Perry County, north of Harrisbuip. in central Pennsylvania and is 
contained within the '^.5-minute Reward quadrangle (fifuie 1 1. The ridpe and vaOey topogr aph y 
of the area, formed from a series of eroded anticlines and syiKiines. is composed of Ordovician 
limestones, sandstones, and shales, running northeast to southwest across the area. Buffalo Moun* 
tain, the ridge examined in this study, runs linearly northeast to southwest, with a maximum height 
of 40*t m. The ndge was chosen for its constant slopes and aspects and uniform forest cover. 
Block diagrams created from the digital terrain data show the topography of the study area ( figures 
1 and .if. Buffalo Mountain has three mafor slope units, moderate im 9" ) and steep (<u 22°) 
rectilinear slcjtes ut northwest aspect and moderate <c. 1 5°» slopes of southeast aspect. The ridge 
is covered hy a mixed deculous oak wiHidland consisting predominantly of white oak. ret* oak. 
chestnut oak. and *'iack oak. A field visit to the study area revealed complete crown closure 
throughout the study area and minimum differences in the cover type, distribution, and density 
between the two predominant aspects. 

Two types of data were integrated to provide the data base used in this study; namely, Landsat 
multisi ectral scanner (MSS) data and digital terrain data. The image-processing system used for 
this study was the IDIMS (Interactive Digital Image Manipulation System) produced by E.S.L. 
(Electromagnetic System Laboratory-) Inc; the particular installation is at NASA/Goddard Space 




Fl«hl C'mter (GSFO. C;rc«nhdt. MD. USA. It » compriied of a Hewlett-Faduird 3000 mini- 
computer with a proframmaMc array pmoemor ami specific software for imapr analysis. The 
various data planes used in the study were created as imape files on the IDIMS system. 

The Landsat Jala chosen lor this study were for July with a solar elevation «»f 55® 

This hiph-sun-angie image was selected hes*ause it depicts the period of maximum foliage, which 
enables optimum discnmination l Williams el al l‘*‘'<*> Furtbermort. it depicts a time close to 
the summer s*»lsiice. when ih-.* topographic effect is at a minimum (Holhen and JustK'e and 

IS mmt likeh . therefore. u> . chosen for analysis of cover types. From the raw Landsat data the 
MSS ' 5 and o 5 band-ratio images were created. Fc' ease of integration the Landsat data were 
resampled to 50 m and legislertd to the I :.4.000 topographic sheet. 

The digital terrain data used in this study were the 30-m Digital Flevation Model (DFM> 
data obtained from the I'.S (Ecological Survey (USCES) Digital .\ppiications Team. Reston. Virginia 
(tigufe 4|. The DFM data should be distinguished from the 200-foot Dc/ense Mapping A^eiwy data, 
also available from the US(tS. The coarser resolulK>n of the latter has been criticised with reference 
to the usefulness of the hon/ontal resolution and the height accuracy (Stow and Fstes I**^*ft. 

The DFM data are derived from orthophotos and are available for selected ''.5-minute quadrangles 
within the l^SA (McFwen and Flassal I ^'' 81 . 

The height accuracy of the DFM data was evaluated by comparing *>1 elevation control points 
on the topographic map with the corresponding height values for the same locations from the 
tenain data before the data resampling. The average deviation of the DEM height values from the 
corresponding topographic map values was 5.92 m. Further examination of the best matching 
point within a 3 X 3 window of the DEM data reduced the average variation to 2.79 m. The 
DEM data were well within the 7-m height accuracy claimed by the USGS (Elassal 1980). From 
the 50-m resampled DEM elevation data, both slope (figure 5a) and aspect (figure 5b) were caF 
cuiated for each elevation point, by use of a four-cell elevation matrix, by applying a modification 



of an al^rithm developed hy Shaipnack and Akin 1 1*>69|. Software was developed to use dope, 
aspect, solar elevation, and a/imuth to calculate the incidence angle for each pixel (figure 5c>. 

The Landsat data were registered to the '.5-minute topographn: map for the Reward Quadrangle 
and the eight surrounding map quadrangles. Fighty-four evenly dispersed control points were 
chosen for the map sheets and the 1 andsat data. The root mean square (rms) error was 0.79 pixeh. 
and tne maximum error was I ."’5 pixels. The five control points located within the Reward Quad- 
rangle had a rms of pixels, with a maximum error of 1 .03 pixek. 

3. Techniques .xaniincd for reducing the topographic effect 

Four techniques for reducing the topographic effect are examined in thb study, namely, 
spectral-hand ratioing and application of a Lambertian model, a modified Lambertian model, and 
a non-Lamhertian model. 

Spectral-hand ratioing in its simplest form consists of dividing the radiance values in one 
channel hy the corresponding radiance salues in a second channel. Ratioing has been proposed 
as a way of reducing mult nlicative effects within multispectral data. e.g.. by Kriegler et aJ. ( 1969). 
Crane t i97| i. and Vincent ( 1 973). The topographic effect of direct light is one such multipli- 
cative effect. A detailed account of hand ratioing as a technique for reducing the topographic 
effect on ground-based radiometer d.ita is presented by Holhcn and Ju.stice ( 1980). who showed 
that ratioing reduced the topographic effect on their data by an average of 8.^ per cent. However, 
a study by Williams cl al. ( I9"’9) on Landsat data of central Pennsylvania showed the spectral- 
band ratioing did not remove slope- and aspect-induced variations to enai'le improved classification. 

The second technique examined was the application of a Lambertian reflectance model. This 
model is described by Justice and Holben (1979) and is based on the assumption that the surface 
being sensed is Lambertian. A Lambertian surface is one that scatters light equally in all directions. 
The radiance from the surface can therefore be modelled by the cosine of the incidence angle, 
where the incidence angle is the angle between the surface normal and the solar beam (Robinson 
1965). 
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Lambertian Model: L*Ln (cot i), 

where L*radiance. 

Ln*nofmalized radiance. 

(■incidence angle. 

The Lambertian model rttumes the reflected light to be directly proportional to the incident 
light. The relationship between data derived with the Lambertian model and radiances receded b> 
ground-based sensors was examined by Holben and JustK'e ( 1970). there was a strong statistical 
relationship bet >een the data sets. The Lambertian model has been examined for correcting sun 
angle and topographic effects in Landsat data (Ckone rt at. 1977. Justice 1978. Hoffer et at. 
1979, Tom era/. 1979. Smith era/. 1 980. Strahler era/. 1979). Hoffer er a/. ( 1 979) applied the 
Lambertian model to Landsat data and found that this correction slightl>' decreased classiflcation 
accuracy for a range of forest cover types. Smith et at. i\ 980) found that the Lambertian assump- 
tion was only valid over a small range of incidenc'e angles for pine-forest cover types. Cicone 
et at. ( 1977) found high statistical correlations oetween Landsat radiances and incidenc'e angle, but 
found that the Lambertial model overcorrected the data. T«> compensate for this overcorrection. 
ChTone et al ( Id**?) mcHlified the Lambertian assumption by multiplying the slope angle by the 
solar zenith angle and then calculating the incidence angle. This modiflcation was determined 
empirically with data for two solar elevations, and it led to improvements over the L *mbertian 
model. 

All natural surfaces have preferred orientations of scattering (Knebel 1976. 1978). and thus 
they are more satisfactorily modelled by a non-Lambertian model. Smith et at. ( 1 980) applied 
a non-Lambertian model to Landsat data, using a function developed by Minnaert (1941. I % I ). 
The model takes into consideration the incidence angle and the exitance angle made with the 
surface, as a function of K. a constant value. 
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Non-lamhertian Model: 


L«Lfi (cos*^i) Icos*^"* e), 
where: L=mdiance, 

Ln- radiance when i-e*0, 

1= incidence. 
e=e\itance. 

K^Minnoen con$Unt 

I he K \j\uc IN ilcrncd fri»rn ihc data anti in eiiui'alcnt to the dope of the reitrejision line formed hy 
pli«tttnic IttfM < L M Ci>N c > I aeaitiNi loe < i ct» 1 1 i con c ^ I The K value gives a measure of the Lambertianess 
of a Niirface. a k ot I MgniliCN a lamK^rtun Nurface. Mod natural surfaces have values betweea 
0 and I Smith er ul \ h>HOi showed that the nt>n*l aniK-rtian model performed significantly better 
ihan the Lambertian nu^del in reducing timin'indiiced variations on pine woodland. Justice and 
HolK-n i examined the Minnaerf tunciion with reference to ground-based data and fcHind a 
similar impnnemenf over the I amhertian rii«Hlel tor reducing the topographic effect. However, sub- 
si'quenf fhetncficjl exanimahon of the k value b> Che authors has shown that the nHXlel is inappro- 
priate ft'r k values ot erealcr Chan lie in cases in whicli there is a high specfual component. 

4 XiuKsis 

I he aiuivsiN of tfu* data is subdivideil inio finir maior sections I he first sivtum describes 
the nuMsurement of !fie lopt>graphic effec t for the I andsat data taken at s^vlar ekrvatton by 
examination o; i .mdoinlv x-le. led lest Mies I hi s<'>, »»nil s<\ iiun presiMifs a detailed analysis of the 
lopographu effect for six simple iransectN lakcn avf%ws the mountain ridge. The third section 
examines die variabihiv »>f ;he I arufvjt radiance vlata for the six transects, after the application of 
speelral-band rafiome ami ilic 1 amberiian .md lu n*l amK*rtian nomul vaiion models to reduce 
the topt>grapinc cMcct The fin.il s<\(um includes appluation of Ihe optimum nonnali/alion pro- 
cedure fi’ die I .tfulN-it vf.iM .iful a comparison o* tiie et^rrecled and fhe uneorreeted radiance data. 



4.1. The topographic effect on selected test sites 

The Landsat data of the Reward Quandrangle for MSS 7 (0.8-1. 1 nm) (figure 6) show the 
topographic effect on Buffalo Mountain ridge, which runs northeast-southwest across the image. 
The deciduous woodland canopy has a high reflectance in the infrared for both dominant aspects, 
but the northwest slope of the ridge facing away from the sun is visibly darker than the southeast 
slope. The MSS 4 and 5 (figure TIfO.S-O.b pen and 0.6-0.7 ptm) show no visible topographic effect, 
because of the low refiectivity of the woodland canopy in these two channels., allowing the 
atmospheric path radiance to dominate. The major variation that can be seen on the mountain 
rii.ge in channels 4 and 5 is caused by sensor banding. This is probably because at low reflectances 
the topographic effect, although present, is obscured by the quantization of the MSS sensor system 
(Holben and Justice. 1^80). The MSS 6 (0. 7-0.8 ^m) image was visually similar to that of MSS 

F or demonstration of the topographic effect on the Landsat data, three type sites were chosen 
to represent the major slope-aspect combinations on Buffalo Mountain, namely, southeast aspect 
and moderate slope, northwest aspect and moderate slope, and northwest aspect and steep slope. 

:st sites were randomly selected for each of the three locations along the mountain ridge. Four 
sites totaling 7(>7 pixels were selected for the southeast aspect location. Four sites totali'-g 660 
pixels and five sites totaling 213 pixels were selected for the northwest aspect moderate- and 
steep-slope sites, respectively. The number of pixels in the latter class was restricted by the extent 
of the slope class. Mean solar incidence angles of 44.5®, 61. 7“. and 73.5“ were calculated for the 
southeast-aspect moderate slopes, northwest-aspect moderate slopes, and northwest-aspect steep 
slopes. C alculation of mean pixel values for each test site revealed distinct grouping in MSS 6 end 7 
with respect to aspect ( figure 8). 

Field checking of the sites before the analysis had shown no observable differences in *hc 
woodland cover for the three locations. The differences in the mean pixel values were therefor 
associated primarily with the topographic variation. The highest mean pixel values asstxiated with 
the oak woodland cover were in MSS 6 for the southeast aspect. However, the highest individual 
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radiances were fo’- MSS which was quantized to 64 levels, half the quantization level of the other 
three channels. The northwest-aspect moderate-slope radiances for MSS 6 and MSS 7 were higher 
than the steep-slope radiances for the same aspect. Three distinct ranges of mean pixel values were 
obtained for the three slope-aspect classes for MSS 6 and 7. For MSS 6 a difference of 1 7 pixel 
values separated the mean values for the sites associated with the southeast and northwest aspects. 
Such a large range indicates that selection of training sites for cover classification from one slope- 
aspect location may not adequately describe the radiances from the same cover from other loca- 
tions. The degree to which this nu y affect the classification results depends largely on ihe location 
and distribution of other cover-class radiances within the classification-feature space. 

4.2 Analysis of sang le transects 

Sjx s;miple transects of L andsat data across the mountain ridge were taken to examine the 
topographic effect in detail and to assess the statistical relationship between the Landsat radiances 
and the data denveii w»ih the pn>posed models, fach transect contained 20 pixels of Landsat 
d;iia I he mean pixel value. sLinujrd deviation, and range were calculated for each transect (table 
1 ). Lxamination of these values \lu>ws the unitormitv oi Landsat radiances K^tween the different 
transects In MSS o. I ranged 3 liaJ the greatest range. 21 pixel values, and in MSS 7. Transect P 
had the gre atest range. 14 pixel values I levation, pixel values, and ineidenee angles for a sample 
iiansect are presented in figure There is a negative rel itionship between pixel value and inci- 
dence angle, the lowest pixel values were associated with the highest incidence angles and vice 
versa. 

The statistical relaliimship between the radiances and the data derived from the models was 
assessed! for each transect vviih Pearsim’s product-moment correlation ctxMTicient (table 2). The 
correlation coefficient gives an mduution o( the stamgth of the linear relationship between two 
variables. ( orrelaiion coeflu ienls tr> were insignificant at the 0.05 level for MSS 4 and 5 for all 
three models Strong positive relationships (r > 0.84) were found for data derived with Ihe Lam- 
bertian and ni>?vl amhertian models, and a strong negative relationship was found for data derived 
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with the modified Lambertian model for all transacts. There was no substantial difference in the 


size of the r values among the three models. 

4.3 Evaluation of topographic eff<'ct reduction techniques 

Spectral-band ratioing and the three models were assessed for reducing the topographic effect 
on Pennsylvania Landsat data. 

Coefficients of variation (CVs) were calculated for the ratios of MSS 6/5 and MSS ~/5 for 
each transect and are presented for comparison with CVs for the raw data (table 3). The CV 
is the standard deviation divided by the mean value, and it is used to compare variations within 
data sets. In this study a decrease in topographic effect is signified by a decrease in the CV. I'hc 
CVs indicated a greater variation in the MSS 7/5 ratio !han in the MSS 6/5 ratio, because of the 
higher variation in MSS 7. The CVs for MSS 6^5 and 7/5 ratios decreased slightly from the CVs 
for the raw MSS data for three out of six transects. The poor reduction in the topographic effect 
by the MSS 6 '5 and 7/5 ratios is due to the negligible variation in MSS 5. The relatively constant 
radiance in MSS 4 and MSS 5 is associated with the high absorption by the green leaves of '.he 
deciduous wocKlIand and means that ratioing is effectively dividing MSS 7 and MSS oy a con- 
stant value. 

Application of the Lambertian model to the data to reduce the topographic effect (i.e.. di- 
vision of each pixel by cos i) led to a large increase in the ('V for each transect over the variation 
for the raw data (table 3>. Dividing the radiance from a Lambertian surface by cos i normalizes 
the radiance to the equivalent radiance from a fiat surface with the sun overhead. This increase in 
the CVs is presumably due to the inapplicability of the Lambertian assumption for the woodland 
cover type in queftion. The largest CV was found for normalized MSS 4 and 5 data, probably 
because of the very small proportion of the MSS 4 and 5 radiance that was directly related to the 
variation in the incidence angle, i.e.. direct light reflectance. Over the majoiity of slopes, the 
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diffuse-light component appears essentially constant (Justice and Holben, 1^80). The modified 
Lambertian model developed b\ ('icone Ji al. ( was a considerable improvement over the 
Lambertiai^ nuxlel m reducini* the tofH>eraphic effect (table 3>, but for all transects, it increased the 
toptvgraphy-indiiced variation ibi>ve that found in the raw Landsat data. 

Application of the non-Lamherliaii model to the transect data involved calculation of the 
Minnaert K .alues for each transect ( table 4). 1 he regression lines used to calculate the K values 
for MSS 4 and 5 had r' values of less than 0.50 and were deemed invalid for any correction of 
the data. The \oy^ r~ values were a direct result of the negligible variation in the MSS 4 and 5 
radiances. The regrv'ssion lines ust‘d to derive the K values for MSS t and ^ had r^ > 0.74, with K 
values ranging from 0.2b to 0 37 tor MSS and 0.20 to 0.32 for MSS (>. The r' values for these 
data w.re higher than tiu>se repc^rted b> Smith i f al { 1^801 and were less variable within the 
cover lypiv 1 he values were then substituted m the non-Lambertian mcK^el |(cos*^i> (cos*^ *e)) for 
each transect, and the l As were calculated ?or the corrected data (table 3f The CVs for the iion- 
L*imbertian moiiel were markcvilv smaller than the C A \ t(>r raw MSS b and data. The maximum 
reductu>n m the C A \ was trom 1 0.:^ to y MSS o and from 1 5.0 to 4.0 tor MSS 

4 4 .\pplic.ition and \N^^.•^smi■nl oi tlu- 'nMi-l amberfu.n im>ilel to correct tlic l,ands,il data 

I !u* tnM’ 1 ,imt'crT...n novle! ,.s slunv?': !o !\- ifu* corievlu^n prov,cdurc for reducing the 
tt^pograpluc ctleci (table » Before .ippK n^g the non-l amlHTtian Muntel fi> the Landsat uala for 
Buffaft> Mountain, it was iu*wessjr\ \o deckle which k value to subsinuie inti> the model. Lse of 
the nunlel with multiple k \alucs would be too complex tor any i>pe (lonal study. Therefore a 
sensitiMtv test was ufulcrtaken assess Inm application i>l a smgle k value to all the transects 
lor MSS <s would ettect the ( A \ allies i iH'tficients oi variation were calculated for each tran- 
sect with data domed liom the non-Lamberiian im>del and K values from 0 to I (table 5). Table 

5 shows considerable vanatuMi m (A lor each transect o\er a range ot K values. A mean k value of 
0 2b'^^ was selected Iriun the transect sUid> and applied to the data. I his mean K value led to a 
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negligible change in the CV value from those derived with individual K values from each transect. 
The non-Lambertian model with a K value of 0.269 was applied to the MSS 6 data to produce a 
normalized image (figure 10). The effectiveness of the non-Lambertian model in reducing the 
topographic effect was determined by comparing the separability of radiance data for the three 
dominant slope-aspect classes before and after normalization. Landsat data for three slope-aspect 
sites were extracted; 1 125 pixels for the northwest moderate slope. 231 pixels for northwest 
steep slope, and 1480 pixels for the southwest moderate slope. The mean values and standard 
deviations are presented in table 6 and show that the difference in pixel values between the three 
sites was considerably reduced by the non-Lambertian model. The difference between the mean 
values for the northwest steep slope and the southeast moderate slope was reduced from 13.9 to 
I.*' <86 per cent) in MSS 6 rnd from 9.7 to 1 .2 (87 per cent) in MSS 7. The v.ariation within the 
sites as measured by the standard deviation increased slightly in rhe corrected data by an average 
of 0.5. The results show that the topographic effect was substantial.y reduced with the non- 
Lambertian model. 

5. Summary of results and conclusion 

It is commonly suggested that Landsat cover-classification studies should be undertaken on 
imageiy obtained at high solar elevation to reduce variations in spectral response because of the 
topographic effect. Results from this study show tha» even high-sun-angle data have an apparent 
topographic effect that may confound cove' classification. Examination of the June l‘J76 Land- 
sat imagery' shewed a marked uifference between the radiances associated with the three dominant 
slope-aspect combinations found in the study area. The topographic effect was evident for MSS 
and MSS 6 with a pixel-value range of 21 among randomly chosen transects across the study ridge. 
No topographic effect was apparent in MSS 4 and 5 because of a combination of high absorption 
in those spectral channels by healthly deciduous vegetation and the Landsat quantization procedure. 

Application of four Piethods for reducing the topographic effect to tiie sample transects across 
the test ridge gave the following results: 



Spcctral-band ratioing (MSS ^ 5, MSS 6-5) slightly decreased the variation of the Landsat 
data for 50 per cent of the transects. Because of the small variation in MSS 5, ratioing could not 
eliminate the direci*light topographic effect, as the denominator of the ratio was virtually constant. 
This reasoning also may explain the p(X>r results obtained by Williams at. ( 1^7d) with band- 
ratioing techniijues. 

When applied to the Lands*it data, the Lambertian mod»d increased the topographic 
effect. This degradation of the radiance is due to the inapplicability of the Lambertian assumption 
to model the bidirectional reflectance characteristics of tiie woixlland surface 

The modifieil Lambertian model developed by ('i.'one cf at. ( decreased the variatmn 
produced by applying the Lamivrti in model but gave even higher variances than those found in 
the raw Landsat data. 

The non-Lamhcrfiaii model developed ^y Smith <7 al. \ I^^HO) markedly decreased r*v8b 
per cent) the variation of the Lamlvit data. and. theret\>re. reduced the topographic effect 

/Mltunigh ihc non I amlx*rtkm model did not cimipletely eliminate the topographic effect, 
it considerably rcdiKcil ihc \ an.ilioii 1 he differeiKi. 1 etvkcen the mean pixel valiiCN asstvialed 
with the two extreme slo| e-.opec( Mtes wjn reduced from to \ per cent) in MSS b and 

to \ .l per cenii m MSS " 1 lie maximum v.inatiim iHMween llie mean corrected pixel 

values from the three slopc-.iNpect conheuralh vas I >> for MSS b jiul \ .l foi MSS ", 

In conclusion, this snidv has demonstrated that high-<|Uj|ily digital lerrain data such as the 
IJSCiS l)l M data can be used via modelling to enhance the utility of multispectral satellite data. 
Digital terrain data can be used to develop and test improved radiative transfer models, w'hich, m 
turn, may lead to improved cover clavsilicatit>n of Landvit data 

Alihougli this study has demonstrated the success of a non-Lambert lan model for reducing 
the topographic effect aiul the potential U>r preprocessing Land liata. further development of 
the existing model is required tor application of tlie technique to more C(Mnplex surface-cover 


conditions 
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Table I 

Mean, standard deviation, and range of Landsat pixel values for six sample transects 

across Buffalo Mountain, PA. 


Transect 

Multispectrai 

scanner 

Mean 

Pixel Value 

Standard 

deviation 

Range 

N 

1 

4 

16.1 

0.64 


20 


5 

12.6 

0.68 




6 

60.0 

4.9 

16 



7 

34.2 

3.2 

9 


2 

4 

16.0 

0.83 

3 

20 


5 

12.6 

0.60 




6 

59.0 

4.8 

15 



7 

33.5 

3.7 

13 


3 

4 

16.2 

0.64 


20 


5 

12.3 

0.73 




6 

60.9 

6.0 

21 



7 

34.5 

3.7 

13 


4 

4 

15.9 

0.72 


20 


5 

12.5 

0.60 

■y 



b 

59.3 

4.6 

15 



1 

33.7 

2.5 

8 


5 

4 

15.8 

0.87 

3 

20 


5 

12.6 





6 

59.5 

4.4 

14 



7 


3.1 

12 


6 

4 

15.9 

0.8? 

3 

20 


5 

12.7 

0.79 

3 



6 

61.6 

6.5 

18 



7 


5.3 

14 
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Table 2 

Pearson's product-moment correlation coefficients (r) for Landsat radiances and 
data derived from the radiance models for six transects. 


Transect 

no. 

Miiltispectral 

scanner 

N 

1 

Model! 

3 

1 

4 

20 

-0.39* 

-0.27* 

0.20* 


5 


0.13* 

-0.08* 

0.0 1* 


6 


0.88 

-0.84 

-0.85 


7 


0.89 

-0.88 

0.88 

n 

4 

20 

-0.07* 

-0.02* 

0.30* 


5 


0.26* 

-0.11* 

-0.11* 


6 


0.90 

-0.93 

0.89 


7 


0.89 

-0.87 

0.87 

3 

4 

20 

0.57* 

-0.57* 

0.40* 


5 


0.54* 

-0.49* 

0.40* 


6 


0.93 

-0.84 

0.91 


7 


0.91 

-0.79 

0.90 

4 

4 

20 

-0.20* 

0.28* 

0.24* 


5 


0.57* 

-0.52* 

0.44* 


6 


0.9(> 

-0.90 

0.96 


7 


0.97 

-0.94 

0.96 

5 

4 

20 

i).24* 

-0.29* 

0.43* 


5 


-0 28* 

-0.28* 

0.08* 


h 


0.S4 

— O.bb 

0.85 


7 


0.8 

-0.78 

0.89 

(> 

4 

20 

0.27* 

-0.32* 

0.24* 


5 


0.27* 

-0.39* 

0.2"* 


<> 


0 94 

-*V94 

0.94 


7 


0.9-' 

-0.95 

0.96 


Model I C'osine i 

2 Moditied cosine i 
} Non-Liinbertian model 
* Insignificant at tlie 0.05 level of probability. 


IS 




Table 3 

C oeiritients of variation for the uncorrected and corrected transect data. 
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Table 4 

Minnaert K values and coe>Ticients of' determinatioii (r*) lor the six sample transecis. 


Transect 

Multispectral 

scanner 

channel 

N 

K 

r* 

1 

4 

:o 

-0.04 

0.07 


5 


0.03 

0.03 




0.26 

0.76 


-» 


0.30 

0.«0 

*> 

I 

:o 

0.01 

0.00 


5 


0 07 

0.16 


6 


0.22 

0.83 


”? 


0.32 

0.82 

3 

4 

20 

0.09 

0.34 


5 


0.13 

0.34 


b 


0.32 

0.85 




0.35 

0.84 

4 

4 

20 

0.01 

0.01 


5 


0.14 

0.50 


b 


0.27 

0.94 




0 27 

0 95 

s 

4 

20 

0.05 

0.10 


5 


-002 

O.Ol 


b 


0.20 

0.74 


7 


0.26 

0.82 


4 

20 

0.05 

0.13 


5 


0.05 

0.10 


b 


0.25 

0.89 


7 


0.3' 

0.94 


:o 





f oeJricicnts of variation (CV) for the non ^ 

a ran,: oS™ '.r"'*’ 



MulfispectraJ 

'.hannd 


6 

7 

6 

7 

6 

7 

6 

7 


— ^ cen t) with indicate d K values 

r\ ^ ^ ~ 



10.3 

11.3 

8.9 

8.8 

7.7 

9.6 

10.8 

15.4 


7.8 

8.8 

5.9 

5.8 

5.1 

6.8 

7.2 
11.8 


5.6 

6.6 

3.2 

3.0 

4.0 
4,7 

4.4 

8.2 


4.3 
5.1 

2.5 

2.3 

5.5 

4.4 

4.5 
5.1 


7.7 

6.3 

11.3 

8.7 

15.2 

11.9 

19.3 

15.5 

23.5 

19.4 

28.0 

23.5 

32.6 

27.2 

4.9 

5.2 

7.0 

6.8 

9.8 

9.2 

12.8 

12.1 

16.0 

15.1 

19.3 

18.3 

22.7 

21.6 

4.9 

4.9 

8.2 

8.1 

11.7 

11.6 

15.3 

15.1 

19.0 

18.8 

27.7 

22.5 

26.6 

26.3 

8.5 

6.5 

12.0 

9.7 

15.7 

13.2 

19.7 

17.0 

23.8 

21.0 

28.1 

25.1 

32.6 

29.4 


7.6 

4.1 


11.6 

6.6 


15.9 

10.4 


20.4 

14.6 


25.0 

19.0 


29.7 

23.5 


34.4 

28.1 


9.8 

10.7 

7.7 
7.5 

7.5 

9.1 

10.5 

15.0 


Tabic 6 

Mean pixel values and standard deviations calculated for uncorrected data and data corrected with the 

non-Lambertian model for three slope-aspect sites. 


Pixel value 


Slope-aspect site 

Mean 

incidence 

angle 

N 


llneorrected data 



Corrected 


MSS 6 

MSS 7 

MSS 6 

MSS 7 

X 

SD 

X 

SD 

X 

SD 


SD 

Northwest (moderate slope) 

6:° 

1125 

56.3 

3.3 

31.9 

2.1 

68.4 

4.1 

40.0 

2.6 

Northwest (steep slope) 

o 

o 

231 

50.7 

2 2 

27.8 

1.3 

67.0 

2.9 

39.7 

1.8 

Southeast (moderate slope) 

72" 

1480 

64.6 

2.6 

37.6 

1.7 

68 9 

2.9 

40.9 

1.8 

Maximum range in means 



139 

9.7 

1.9 

1 . 
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